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New neutron diffraction and inelastic neutron scatter- 
ing experiments on Ba2CuGe207 suggest that the previ- 
ously suggested model for the magnetism of this material 
(an ideal sinusoidal spin spiral, stabilized by isotropic ex- 
change and Dzyaloshinskii-Moriya interactions) needs to be 
refined. Both new and previously published experimental re- 
sults can be quantitatively explained by taking into account 
the Shekhtman-Entin-Wohlman-Aharony (SEA) term, a spe- 
cial anisotropy term that was predicted to always accompany 
Dzyaloshinskii-Moriya interactions in insulators. Our exper- 
imental results present the first clear evidence to that SEA 
interactions can lead to substantial observable effects in a real 
magnetic system. 
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Among the more exotic magnetic interactions in solids 
is the so-called asymmetric exchange, first predicted the- 
oretically by Dzyaloshinskii Q. Unlike conventional 
Heisenberg exchange coupling, that is proportional to 
the scalar product Si ■ S2 of interacting spins, asymmet- 
ric exchange is proportional to the corresponding vector 
product. In the spin Hamiltonian it is usually written as 
D{Si X S2), where D is the Dzyaloshinskii vector associ- 
ated with the bond between the two interacting magnetic 
ions. A microscopic model for asymmentric exchange in- 
teractions was first proposed by Moriya Q , and is essen- 
tially an extension of the Anderson superexchange mech- 
anism 1^, that allows for spin-flip hopping of electrons. 
While forbidden by symmetry in centrosymmetric crys- 
tal structures, Dzyaloshinskii-Moriya (DM) interactions 
were found to be active in a number of non-centric com- 
pounds, where they lead to either a weak ferromagnetic- 
or helimagnetic- distortion of the coUinear magnetic state 
The inclusion of the DM term breaks 0(3) invari- 
ance of the originally isotropic Heisenberg spin Hamilto- 
nian, reducing the symmetry to 0(2): to take full advan- 
tage of the cross product term the interacting spins must 
be perpendicular to the vector D. DM interactions thus 
play a role of an effective two-ion easy-plane anisotropy, 
with the easy plane normal to the vector D. 

Only recently Schekhtman, Entin-Wohlman and 
Aharony realized that there is more to Moriya's mech- 



anism than just the vector-product term For very 

fundamental reasons the DM cross-product must always 
be accompanied by a two-ion easy axis anisotropy term 
that exactly compensates the easy-plane effect of the vec- 
tor product. The additional (SEA) term can to a good 
approximation be written as ■^{SiD){S2D), where J is 
the Heisenberg (isotropic) component of superexchange 
coupling. Often referred to as "hidden symmetry", the 
SEA term restores 0(3) invariance of the Hamiltonian, 
at least locally. Originally, the SEA term was invoked 
to explain the spin anisotropy in the orthorhombic phase 
of La2Cu04 It was later realized that this term 

alone can not account for all the observed effects, par- 
ticularly for the magnetic anisotropy seen in the tetrag- 
onal phase [ p^pHJl^ . To our knowledge, to date there 
has been no "clean" experimental evidence unambigu- 
ously pointing to the presence of SEA interactions. In 
the present paper we present such experimental data for 
the helimagnetic insulator Ba2CuGe207. We demon- 
strate that only by taking into account the SEA term one 
can obtain qualitatively and quantitatively correct pre- 
dictions for the magnetic structure and spin wave spec- 
trum. 

As was shown in a series of recent publications 
]l|,0jl5Ul|],Ba2CuGe2O7 is a particularly useful model 
system for studying DM interactions. The magnetism of 
this compound is due to Cu^+ ions that form a square 
lattice in the (a, 6) tetragonal plane of the crystal. The 
principal axes of this square lattice, hereafter referred 
to as the x and y axis, run along the [110] and [110] 
crystallographic directions, respectively. To complete 
the coordinate system we shall choose the z along the 
[GDI] direction. In the magnetically ordered phase (below 
Ttv — 3.2 K) all spins lie in the (110) plane. The magnetic 
propagation vector is (1 -I- C, Ci 0)i where C = 0.0273, and 
(1,0,0) is the antiferromagnetic zone-center. The mag- 
netic structure is a distortion of a Neel spin arrangement: 
a translation along (i,i,0) induces a spin rotation by 
an angle <j) — 2tt( « 9.8° (relative to an exact antiparal- 
lel alignment) in the (1,-1,0) plane. Along the [110] 
direction nearest-neighbor spins are perfectly antipar- 
allel. Spins in adjacent Cu-planes are aligned parallel 
to each other. Only nearest-neighbor in-plane isotropic 
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superexchange antiferromagnetic interactions are impor- 
tant (J « 0.96 meV, as determined by the measured 
spin wave bandwith The hehcal state is stabiUzed 

by DM interactions. In the current model Dzyaloshinskii 
vectors for nearest-neighbor Cu-Cu pairs he in the (x, y) 
plane and are oriented perpendicular to their correspond- 
ing bonds: D\\y for an cc-bond and D\\x for y-bonds, 
respectively (Fig. 1 in Ref. |^). The corresponding en- 
ergy scale is £> « 0.17 meV. In the discussion below we 
shall use the numerical values for J and C, quoted above 
as given, and perform all calculations without using any 
adjustable parameters. 

The magnetic structure of Ba2CuGe207 is expected 
to be particularly sensitive to the presence of SEA in- 
teractions. Indeed, in Ba2CuGe207 SEA easy axes that 
correspond to the y- Cu-Cu bonds are parallel to the x 
axis, i.e., are in the plane of spin rotation. Regions of 
the slowly rotating spin spiral where the local staggered 
magnetization I is almost parallel to x become more en- 
ergetically favorable than those where I is almost paral- 
lel to z (crystallographic c-axis). SEA anisotropy must 
therefore lead to a distortion of the ideal sinusoidal spi- 
ral, and modify the period of the structure. The SEA 
term is expected to produce exactly the same distortion 
as a magnetic field H applied along the z axis: the latter 
also has the effect of forcing the local staggered mag- 
netization into the {x, y) plane. The role of a 2-axis 
field is rather dramatic and has been studied in detail 
| p5| , p^ . The period of the spiral increases with increasing 
H and diverges at He ~ 2.15 T, resulting in a commen- 
surate spin-flop antiferromagnetic state aX H > He- For 
< H < He the spin structure is described as a "soliton 
lattice" where regions of the commensurate phase are in- 
terrupted by regularly spaced antiferromagnetic domain 
walls. In the soliton phase, in addition to the princi- 
pal magnetic Bragg peaks at (1 ± C, ±Ci O); characteristic 
of an ideal spiral, one expects to see all odd magnetic 
Bragg harmonics at (1 ± 3C, ±3C, 0), (1 ± 5^, ±5C, 0), etc. 
By comparing the experimental field dependencies of C 
and the higher-order Bragg peaks to theoretical predic- 
tions for the "DM-only" (Ref. |l|,|lll) and "DM+SEA" 
models we can hope to obtain direct evidence for SEA 
interactions in Ba2CuGe207. 

We can make the above discussion quantitative by in- 
cluding the SEA term into the phenomenological energy 
functional that was previously used to describe the be- 
haviour of Ba2CuGe207 in the framework of the DM-only 
model ||l^,^. This procedure is rather straightforward 
and the principal conclusion is that all previously ob- 
tained DM-only results can be recycled, by replacing H 
in all equations by the effective field 



H'^^^Jm + 2Aps/{x±-X\\)- 



(1) 



Here ps ~ JS^ is the spin stiffness, x± ^^'^ X\\ ^re the 
local transverse and logitudinal susceptibilities, respec- 



tively, and the SEA term is represented hy A = o? 12 k, 
D'^/2,P. The parameter a is defined by tana = D/J, 
and is equal to the spin rotation angle (j) in the DM-only 
model. According to our continuous-limit calculations, 
in the DM-I-SEA model a = m:cia.\\{D/J) = 

First, let us consider the field dependence of the incom- 
mensurability parameter ^ that for the DM-I-SEA model 
can be obtained by replacing H by H°^ in Eqs. 4 and 7 
in Ref. |l|]: 
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Here (3 is an implicit variable and K and E are complete 
elliptic integrals of the first and second kind, respectively. 
In Fig. 0we replot the C(^) data from Ref. [|l^ in reduced 
coordinates. The solid and dashed lines are the theoreti- 
cal curves plotted with and without taking into account 
SEA interactions, respectively. We see that the inclusion 
of the SEA term hardly affects the shape of the C,{H) 
curve. However, the theoretical prediction for H^. is sub- 
stantially different in the DM-only and DM-I-SEA mod- 
els. Combining Eq. |l] from above with Eq 5 in Ref. |l^] 
one readily obtains: 



Hr 



(4) 



For the low-temperature limit in Ba2CuGe207 we can use 
the classical expressions ps — JS'^ — 0.24 meV, X|| — 
and x± = (gcf^B)^ /8J, where gc — 2.47 is the c-axis gy- 
romagnetic ratio for Cu^+ in Ba2CuGe207 ||l^. One can 
expect these classical estimates to be rather accurate, as 
they rely on the effective exchange constant J, that itself 
was determined from fitting the classical spin wave dis- 
persion relations to inelastic neutron scattering data [ p^ . 
Substituting a = 27r||C = 0.177, we get the estimate for 
the critical field He — 2.05 T. This value is much closer 
to the experimental value He ~ 2.15 T, than our previ- 
ous estimate He ~ 2.6 T |l^], obtained without taking 
into account the SEA term. 

As mentioned, SEA interactions have a substantial in- 
ffuence on the intensity of higher-order Bragg harmon- 
ics. In the DM-only model in zero field higher-order 
Bragg reffcctions are totally absent. For the DM+SEA 
model, combining our expression for H'^^ with Eqs. 17,18 
in Ref. [|l^, for the relative intensities of the 1st and 3rd 
harmonics, in the small field limit (weakly distorted spi- 
ral) \(p — a\ (p we get: 
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In zero field this gives h/h = 1/256 « 4 x 10"^. To 
verify the relation (^ we performed new magnetic neu- 
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tron scattering experiments on Ba2CuGe207 single crys- 
tal samples. The measurements were done in two ex- 
perimental runs, on the IN-14 3-axis spectrometer at 
the Institut Laue Langevein (ILL) in Grenoble, and the 
SPINS spectrometer at the Cold Neutron Research Fa- 
cility at the National Institute of Standards and Tech- 
nology (NIST). The samples were similar to those used 
in previous studies In each experiment the crystals 
were mounted with their c-axes vertical, making (h, k, 0) 
wave vectors accessible for measurements. The magnetic 
field was produced by standard split-coil superconducting 
magnets. The sample environment was either a pumped- 
^He or pumped-^He cryostat. The data were collected at 
temperatures in the range range 0.35-5 K. Neutrons of 
energies 3.5 meV or 2.5 meV were used in most cases. A 
40' — — 40' — A — 40' collimation setup was used with a 
Be filter positioned in front of the sample. In Fig. §(a,b) 
we show some typical elastic scans along the (1 + e, e, 0) 
reciprocal-space line measured in Ba2CuGe207 at low 
temperatures in zero and H = 1 T applied fields. Even 
in the zero-field data in addition to the Ist-order prin- 
cipal magnetic refiection one clearly sees the 3rd order 
harmonic. The measured field dependence of /3//1 (ra- 
tio of Q-integrated intensities) is shown in Fig. ^(c) . In 
our measurements we have taken special care to verify 
that the relative intensities of the two peaks are totally 
independent of the T—H history of the sample. The solid 
and dashed lines in Fig. ||(c) represent the predictions of 
the DM-FSEA (Eq. |) and DM-only (Ref. mod- 
els, respectively. For these theoretical curves we used the 
same numerical values as in our estimates for iJc, and 
no adjustable parameters. An almost perfect agreement 
between the DM-I-SEA model and the experimental data 
is apparent, and so is the failure of the DM-only model. 

It is clear that the SEA anisotropy term, in addition 
to modifying the ground state, will affect the spin wave 
spectrum. For an ideal spin spiral (DM-only model) the 
classical spin wave dispersion relations can easily be ob- 
tained analytically using the Holstein-Primakov formal- 
ism. The result of such a calculation for Ba2CuGe207 
is shown in Fig. ||(a). Two acoustic branches (hereafter 
referred to as the ±C modes) emerge from the two mag- 
netic Bragg peaks at (1 ± C, CjO)- ^ third branch (the 
0-mode) has a gap at the antiferromagnetic zone-center, 
equal to the Dzyaloshinskii parameter D. This branch al- 
most excatly passes through the intersection point of the 
±C modes. The actual dispersion curves in Ba2CuGe207 
were measured in constant- Q inelastic scans using the ex- 
perimental setups described above in the fixed-incident- 
energy mode. In these experiments the typical energy 
resolution, as measured by scanning through incoherent 
scattering, was 0.075 meV FWHM at zero energy trans- 
fer. A strong incoherent signal and the "tails" of mag- 
netic bragg reflections prevented us from collecting reli- 
able data for energy transfers of less than w 0.17 meV. 
A typical inelastic scan (raw data) is shown in Fig. ^. 



Combined data from the two series of experiments are 
summarized in the experimental dispersion relations in 
Fig. ||(b) (symbols). Two prominent features of the mea- 
sured dispersion curves are to be discussed here, i) The 
two ±C modes do not intersect at the Neel point. In- 
stead, at Q = (1, 0, 0) there is a clear repulsion between 
these two branches. This repulsion is again manifest at 
Q = (1 -I- 2C, 2C, 0) and is seen as a discontinuity in the 
+C- branch, ii) At the antiferromagnetic zone-center the 
0-mode lies visibly lower than the extrapolated point of 
intersection of the ±C branches (dashed lines). Obvi- 
ously, the DM-only model fails to reproduce the observed 
dispersion relations. An exact analytical calculation of 
the spin wave spectrum in the presence of SEA interac- 
tions is not possible. However, in the limit \<f> — a\ 4>^ 
a condition well satisfied in Ba2CuGe207 in zero applied 
field, one can use the series expansion method described 
in Refs. 20|, and leaving only the lowest-order terms 
in A. Performing this somewhat tedious calculation for 
Ba2CuGe207, we obtain the dispersion relations shown 
in solid lines in Fig. ^(b), and find very good agreement 
with experiment with no adjustable parameters. 

In summary, both the static and dynamic magnetic 
properties of Ba2CuGe207 contain quantitative evidence 
of the SEA anisotropy term. The issue that we tried 
to address above is not whether or not SEA interactions 
exist: if one believe Anderson's and Moriya's superex- 
change mechanisms, one is forced to accept the presence 
of the SEA term as well. Rather, our results for the first 
time show that SEA interactions can result in very in- 
teresting measurable effects in a real magnetic system. 
We also would like to emphasize that the new results 
do not undermine our previous conclusions regarding 
the commensurate-incommensurate Dzyaloshinskii tran- 
sition in Ba2CuGe207. The transition occurs exactly as 
described, only now we have a more accurate theoretical 
estimate for the critical field He- 
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FIG. 4. Typical inelastic scans measured in Ba2CuGe207 
at T = 1.5 K and T = 3 K (combined data in lower plot). 
The shaded curves represent the individual Gaussians in a 
multi-peak fit (heavy solid line). The gray area shows the 
position of a "Bragg-tail" spurious peak. 



FIG. 1. Field dependence of the incommensurability pa- 
rameter (, as previously measured in Ba2CuGe207 (Refs. 
[ p"5|Jl^ ). The solid and dashed lines are theoretical predic- 
tions that do (this work) or do not (Ref. jl^) take into ac- 
count SEA interactions, respectively. 



FIG. 2. Typical elastic scans along the (1,1,0) direction 
in the vicinity of the antiferromagnetic zone-center (1,0,0), 
measured in Ba2CuGe207 at T = 0.35 K in zero field (a) and 
in a = 1.6 T magnetic field applied along the c-axis (b). 
(c)The square root of the measured ratio of the intensities of 
the (1 -|- 3^,3(^,0) and (1 + CiCjO) peaks plotted against the 
square of the applied field. The lines are guides for the eye in 
(a) and (b) and theoretical curves in (c), as in Fig. ^ 



FIG. 3. (a) Classical spin wave dispersion relations cal- 
culated for Ba2CuGe2 07 without taking into account SEA 
interactions, (b) Solid lines: same, with the SEA term in- 
cluded. Symbols: experimental dispersion curves measure in 
Ba2CuGe207 at T = 0.35 K and T = 1.5 K with inelastic 
neutron scattering. 
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Zheludev et al. Fig. 1 
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